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Abstract 

In this work we consider triangulations of point sets in the Euclidean plane, i.e., max- 
imal straight-line crossing-free graphs on a finite set of points. Given a triangulation of 
a point set, an edge flip is the operation of removing one edge and adding another one, 
such that the resulting graph is again a triangulation. Flips are a major way of locally 
transforming triangular meshes. We show that, given a point set S in the Euclidean plane 
and two triangulations T± and Ti of S, it is an APX-hard problem to minimize the number 
of edge flips to transform T\ to T^- 

1 Introduction 

Given a finite set S of n points in the Euclidean plane, a triangulation T of S is a maximal 
straight-line crossing-free graph on S. An edge flip is the operation of removing an edge e 
of T and adding a different edge / such that the resulting graph T is again a triangulation of 
S. This requires the two empty triangles incident to e to form a convex quadrilateral, which 
is the same as the one formed by the triangles at / in T. The flip operation defines the graph 
of triangulations Q of S, also called the flip graph of S. For a given set S, the vertex set 
of Q is the set of all triangulations of S. Two vertices in Q are adjacent if the corresponding 
triangulations can be transformed into each other by a single edge flip. Lawson [T7] showed 
that Q is connected for any S with diameter 0(n 2 ). Hurtado, Noy, and Urrutia [15] proved 
that this bound is tight. 

Bose and Hurtado [8] give an extensive survey on the flip operation. Flips in triangulations 
are used for enumeration and as a local operation to generate meshes of good quality according 
to a predefined criterion. For example, Lawson [18] showed that one can always obtain 
the Delaunay triangulation after 0(n 2 ) locally improving flips. The Delaunay triangulation 
optimizes several criteria. Also, heuristic methods for improving other properties of triangular 
meshes may apply local optimization using flips in combination with techniques like simulated 
annealing. See [2] for information on the topic of mesh optimization. Another reason for 
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the continuing interest in flips in triangulations is the bijection between binary trees and 
triangulations of convex point sets. There, a flip corresponds to a rotation in the binary 
tree. Properties of the flip graph for convex point sets were studied in the landmark paper of 
Sleator, Tarjan, and Thurston |22j . 

For general point sets, Hanke, Ottmann, and Schuierer |13j show that the length of a 
shortest path between two triangulations in Q (i.e., the flip distance) can be bounded from 
above by the number of crossings between the edges of the two triangulations. Eppstein [12] 
gives a polynomial-time algorithm for computing a lower bound; note that the point sets for 
which Eppstein's result is tight may not contain empty convex 5-gons, which requires more 
than two points being placed on a common line for sets of 10 or more points (see, e.g., [1]). 
Throughout this paper, we make the common assumption that S is in general position, i.e., 
that no three points are collinear. 

Despite these results, the complexity of determining the flip distance between two trian- 
gulations has been open. Our main result is that the problem is APX-hard, which sheds light 
on a "fundamental open issue" [8] in the study of flip graphs. It has been addressed as an 
open problem in |13| already in 1996, and, most recently, in a monograph by Devadoss and 
O'Rourke p. 71]. APX-hardness of the problem implies that no polynomial time approx- 
imation scheme (PTAS) exists (i.e., there is no polynomial time algorithm that approximates 
the flip distance by a ratio of at most 1 + e for every constant e > 0), unless P = NP. Most 
recently, NP-completeness of the problem has simultaneously and independently been shown 
by Lubiw and Pathak (19]. However, their reduction is from the Planar Cubic Vertex 
Cover problem, for which a PTAS exists [5j [6] (see also [H p. 369]), and the reduction 
can therefore not be adapted directly to show APX-hardness. For triangulations of simple 
polygons, Aichholzer, Mulzer, and Pilz |2] recently showed that the corresponding problem 
is NP-complete. Interestingly, the question is still open for point sets in convex position 
(or equivalently, convex polygons), regardless of the intensive investigation of that structure 
within the last 25 years (see, e.g., [201 [21] f° r two exemplary results, one long-standing, the 



other recent). A previous preprint version of this paper (arXiv: 1206.3179vl) only showed 
NP-completeness of the corresponding decision problem. 

Clearly, the problem is an NP optimization problem. Our reduction is from the well-known 
Minimum Vertex Cover problem on cubic (3-regular) graphs. In the next section, we show 
certain properties of triangulations of a class of point sets called double chains, which will 
be subsets in our construction. In Section [3[ we present the gadgets used in our reduction 
and analyze the construction. In that section, we only present a rough overview on how the 
points of the set are placed, a more detailed description of how to calculate their coordinates 
is given in Appendix [A} 

2 Double-Chain Constructions 

A main ingredient of our reduction will be gadgets consisting of subsets that, being considered 
on their own, would require a number of flips that is quadratic in their size. 
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Figure 1: A double chain. The points are divided in an upper and lower chain, each chain 
being in convex position in a way that every point of the lower chain "sees" every vertex of 
the convex hull of the upper chain, and vice-versa. 

2.1 A Single Double Chain 

The construction shown in [15] to have quadratic flip graph diameter is the so-called double 
chain. In a double chain D, there are 2n points, n on the upper chain and n on the lower 
chain. Let these points be (u\, . . . , u n ) and (li, . . . , l n ), respectively, ordered from left to right. 
The upper chain is reflex w.r.t. any point of the lower chain, and vice- versa. See Figure[TJ Let 
Pd be the polygon (l\, . . . , l n , u n , . . . , ui). The edges lijUj+i and kk+i for 1 < i < n have to 
be part of every triangulation of D, since there does not exist a straight-line segment between 
two points of D that crosses any of them (such edges are called unavoidable). Therefore, we 
only need to consider the triangulation inside Pd for the following result. 

Theorem 1 ([E]). Consider any triangulation T\ of D where u\ is adjacent to each of 
h, . . . , l n , and any other one, called T2, where l\ is adjacent to iti, . . . , u n . The flip distance 
between T\ and T2 is at least (n — l) 2 . 

See Figure [2] for the relevant parts of the two triangulations. In their proof, Hurtado et 
al. [T3] label the triangles that have two points on the upper chain with 1 and the ones with 
two points on the lower chain with 0. Consider a horizontal line £ that separates the two 
chains. The triangles crossed by £ from left to right define a sequence a of (n — 1) elements 
labeled and (n — 1) elements labeled 1, see Figure [3j Note that there are no triangles 
of a third type stabbed by £. Further note that we do not care about the triangulation of 
the convex hull of either chain; the lower bound on the flip distance stems from the part 
stabbed by £. It is easy to see that only an edge with two differently labeled triangles can 
be flipped in the stabbed part. This corresponds to exchanging an adjacent pair of and 1. 
Flipping the first triangulation to the second one corresponds to transforming the sequence 
0-1 = ((O)™- 1 ^)™- 1 ) to a 2 = ((l) n - 1 (0) n - 1 ), which leads to the desired bound. 

Our next step will be to gain more insight into the way the flip graph is altered by the 
addition of points. For the following definition refer to Figure [4} 

Definition 1. Let D be a double chain of n points, and let W\ be the double wedge defined 
by the supporting lines of u\U2 and l\l 2 whose interior does not contain a point of D. W n 
is defined analogously by the supporting lines of u n u n -\ and l n l n -i- Let W = W\ U W n be 
called the wedge of the double chain. The unbounded set W U Pd therefore is defined by four 
rays and the two chains. A point is outside of D if it is not contained in W U Pd ■ Let £ be 
a horizontal line that separates the two chains. All points in the same half-plane of £ as u\ 
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Figure 2: Two (partial) triangulations of the double chain with a flip distance of at 
least (n - l) 2 . 




Figure 3: An illustration of the labeling argument for the lower bound. By the flip, the 
sequence changes from (11000101) to (11001001). 



are called the upper points. The lower points are defined analogously. The kernel of a double 
chain D is the intersection of the closed half-planes below u\Ui and u n -\u n , as well as above 
hh and l n -il n Y\ 




Figure 4: The polygon and the wedge (gray) of a double chain. The diamond-shaped kernel 
can be stretched arbitrarily by flattening the bend of the chains. 

Let us add a point s inside the kernel of D. From any triangulation of the resulting 
set D U {s}, we can flip the edges between the chains such that they are incident to s. 
Reaching this canonical triangulation only requires a linear number of flips. This behavior 
is well-known, see, e.g., [23J. Consider the case where s is placed outside of Pr> but inside 
the kernel of D. Add edges from s to u± and l\ to again have a triangulation, as shown in 
Figure [5j Then, for flipping all possible edges to be incident to s, we need at most 2n — 2 
flips. 

Proposition 2. Let S be any set of points outside of D. Then, the diameter of the flip graph 
of S U D remains in f2(n 2 ), where \D\ = 2n. 

1 Note that the kernel of D might not be completely inside the polygon Pd (but no point in the kernel is 
outside D). This is in contrast to the common use of the term "kernel" in visibility problems for polygons. 
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Figure 5: An extra point in the kernel of D allows flipping one triangulation of Pd to the other 
in An — 4 flips. Note that an edge common to source and target triangulation is temporarily 
flipped. 




Figure 6: Mapping a triangulation to a local triangulation of a double chain. To the left, 
all triangles intersecting with W are shown. Visually, one can think of "cutting" the edges 
at the boundary of Pd U W (middle) and moving (and merging) the endpoints to the next 
point (right). 

In order to prove the proposition, we consider a mapping L from the set of triangulations of 
S U D to the set of triangulations of the polygon Pd . When flipping an edge in a triangulation 
T of D U S, at most one edge is flipped in the corresponding triangulation L(T) of Pd- 

Consider any triangulation T of S U D. If all edges of Pd are present, L(T) equals the 
triangulation of Pd in T (note that this is also the case for the triangulations T\ and T2 of 
Theorem [I]). Otherwise, consider the following construction (see Figure [6] for an example). 
For any edge e of T that intersects W, we draw a new edge e' of L(T) in the following way. 
If one of the endpoints of e is on a vertex of Pd, e' has the same endpoint. If e passes 
through an edge tijtij+i or ljlj + i, then the corresponding upper or lower endpoint of e' is set 
to itj.fi or lj+i, respectively. If e passes through one of the rays defining W U Pd, e.g., the 
one through u n , then the upper endpoint of e' is placed at u n , such that e' is contained in 
Pd- Let T' = L(T) be the graph induced by the new edges, and let the edges of T that pass 
through W but do not have an endpoint in D be called wide edges. We call the construction 
T" the local triangulation of D when T is clear from the context. The following lemmata show 
that T 1 actually is a triangulation. Observe that, informally, the construction corresponds to 
continuously introducing the edges of the chains along the arrows drawn in Figure |6| while 
continuously sliding the edges of T accordingly. 
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Lemma 3. For every wide edge e G T that is mapped to e' G T', there is an edge e £ T that 
has an endpoint p £ D that is also mapped to e' . 

Proof. Let e' be Uilj. Consider first the case where both endpoints of e' are on the same side 
of (the directed line supporting) e. Consider the empty triangle t of T incident to e that has 
its apex a on the same side of e as e' . If a is outside of D, then another wide edge / of t is 
also mapped to e' . In that case we continue the argumentation with /, as e and / are both 
mapped to the same edge. If a is not outside, then a equals either u% or lj, as otherwise t 
would contain one of them. Hence, one of the edges of t incident to a is also mapped to e' . 

For the case where the two endpoints of e' are on different sides of e (i.e., one of the 
endpoints of e' is u n or l n ), the argumentation is almost the same. W.l.o.g. let i = n and lj 
be to the right of e (note that j might be n). Therefore, u n is to the left of e. Again, consider 
the empty triangle t of T incident to e with apex a to the right of e. Again, if a is outside of 
D, there is another wide edge / of t that is also mapped to e'. If a is not outside of D, then 
a = lj; this follows from the construction of D and the fact that the lower endpoint of e is 
outside D. Hence, an edge of t incident to a is also mapped to e' . □ 

Lemma 4. Every point p of D is incident to at least one edge e ofT such that e disconnects W . 

Proof. This follows directly from the construction of D: If there were no such edge, then there 
would be an angle larger than ir incident to p, and the wedge defined by this angle would 
contain points. This would contradict the fact that T is a triangulation. □ 

Lemma 5. T' is a triangulation of Pp. 

Proof. We have to prove that T' is crossing- free and maximal in Pjj. 

Lemma [3] allows us to only consider non-wide edges. With all relevant remaining edges 
of T being incident to a point in D, the fact that T' is crossing- free follows from T being 
crossing- free, as the mapping only "moves" the endpoints of the edges of T to the next point 
of D. 

If T 1 were not maximal, there would exist a quadrilateral q inside Pjj that is spanned by 
points of D and whose interior does not intersect any edge. If q is not convex, this would 
mean that no edge of T is incident to the reflex vertex of the quadrilateral. But this can not 
happen due to Lemma [4j If q is convex, it is of the form lili+iUj+iUj . If an edge of T would 
have passed through the side kuj, the quadrilateral would not be empty of edges. Hence, kuj 
must have been a part of T. Also, neither liUj+i nor Ujli+\ was present in T, and therefore, 
there was at least one edge crossing the segment /jiij+i, as well as one crossing Ujk + i. None 
of these edges was incident to U or Uj, since otherwise q would not be empty. If there have 
been edges that have crossed both these diagonals of q, choose the leftmost one. As in the 
proof of Lemma [3j this is a wide edge that is adjacent to a triangle with either Zj or Uj. One 
of the other triangle edges would be mapped to kuj + ± or Ujk + i, respectively. But if there 
were two different "leftmost" edges that prevented visibility along the diagonals, these edges 
would have crossed in T (since none of them has crossed l{Uj or was incident to one of huj), 
which contradicts the fact that T was a triangulation. Hence, there is no empty quadrilateral 
in Pd, which completes the proof. □ 
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Since any flip in T results in at most one edge being flipped in T", the lower bound 
construction holds: a shorter flip sequence with outside points would immediately imply a 
shorter flip sequence between T\ and T2 in the proof of Theorem [T] This completes the proof 
of Proposition [2] 

2.2 Multiple Double Chains 

Proposition [2] is, however, of little use when we try to construct a point set that contains 
many double chains and try to argue that the flip distance between two triangulations of the 
set is bounded by the sum of the distances between the local triangulations of these double 
chains. One could imagine that a flip in the overall triangulation leads to changes in the local 
triangulations of several double chains. In this section, we prove that this is not possible. 
Keep in mind that it is a necessary condition that, for any double chain, all other double 
chains are outside and their polygons do not intersect. 

Lemma 6. Let e be a flippable edge of any triangulation T of DuS that is mapped to the edge 
e' in the corresponding local triangulation T' of a double chain D. Then flipping e changes 
the local triangulation only if there is no other edge mapped to e' . 

Proof. Suppose e would not be the only edge mapped to e 1 . If we remove e from T, the graph 
on D defined by the mapping is still the local triangulation T' . If we add the new edge / 
after the removal of e, / must also be mapped to some existing edge /' in T' , as otherwise T' 
would not be a triangulation. □ 

Note that because of Lemma [HJ flipping an edge that is wide for a double chain does not 
change the local triangulation of that double chain. Therefore, a flip can only represent a flip 
in at most four local triangulations (which would already be enough for our reduction, as will 
become clear). However, we can actually prove the following, more accurate result. 

Lemma 7. Let D\ and D2 be two double chains in a point set S. If each of D\ and D2 is 
outside of the other and Pjj 1 n Prj 2 = 0j each flip in a triangulation of S affects at most one 
of the two local triangulations. 

Proof. If the flipped edge e or its replacement / would not both have an endpoint in the same 
double chain D, at least one of e or / would either not dissect the corresponding wedge W 
or be a wide edge of D. It follows from Lemma [6] that such a flip does not influence the local 
triangulation of D. Hence, in the only remaining case there is a quadrilateral that has two 
adjacent points in D\ and two adjacent points in D2 and contains a flippable edge. Let the 
quadrilateral be abed. W.l.o.g., let a and b be part of D\ and e = ac. Suppose, for the sake 
of contradiction, that we flip the edge ac and the flip changes the local triangulation of D\. 
Then ac has to dissect the wedge of D\. However, then also ad dissects the wedge and crosses 
the same edge of Pd 1 as ac. Hence, ac and ad are mapped to the same edge in the local 
triangulation of D±, a contradiction due to Lemma [6| □ 

Corollary 8. If a point set consists of m double chains, each of size 2n, and for every double 
chain all other points are outside, then the flip graph diameter of the whole set is in Q(mn 2 ). 
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Figure 7: An embedding of a cubic graph with the (almost straight) circular arcs at each edge 
ending at a fixed distance around each vertex. 

3 The Reduction 

Up to now we have gathered enough knowledge about double chains as sub-configurations 
in order to use them as the main building blocks in a reduction. We reduce from Cubic 
Minimum Vertex Cover, which is known to be APX-complete [3]. 

Problem 1 (Cubic Minimum Vertex Cover). Given a simple cubic (3-regular) graph 
G = (V,E) with n = \V\, choose a set C C V such that every edge in E has at least one 
vertex in C and such that \C\ is minimized. 

We follow the common approach of embedding the graph G and transforming its elements 
to geometric gadgets. The gadgets consist of points together with the corresponding edges 
in the source triangulation T\ and in the target triangulation T2. We give the overall idea of 
how to embed the gadgets, for a detailed description on how to exactly place the points with 
rational coordinates having a representation bounded by a polynomial in the input size using 
polynomial time can be found in Appendix [Aj 

3.1 Gadgets 

Given an instance (G = (V,E),k) of the Cubic Minimum Vertex Cover problem, with 
n = \V\ and m = \E\ (note that m = ^), we place the elements of V as the vertices of a 
convex n-gon and draw the straight-line edges between them (as a mean of orientation for 
our construction of the point set). Let x' be the maximum number of crossings a single edge 
has, and let x = [y] . It is easy to see that x' < m — 5. For each edge e mark a point c e E e 
that has no more than x crossings on either side along e. Let t be a vector perpendicular 
to e of sufficiently small length (which will be specified in Appendix |A|). Make two copies of 
e and translate them by t and —t, respectively, to obtain the tunnel of the edge, i.e., the 
quadrilateral defined by the two copies of e. Then slightly "bend" the copies towards the 
midpoint of e to obtain two circular arcs A e and A' e . The endpoints of e have to see any point 
on A e and A' e . See Figure [7J 

3.1.1 Edge Centers 

Cf. Figure [8j Instances of the double chain are the main ingredient in our reduction. Let e be 
a straight-line edge of G, drawn between the points v and v' . In a close neighborhood of c e , 
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Figure 8: The double chain at the center of an edge with the source and the target triangu- 
lation. 

place a double chain D e , the edge center, of 2d points (we will fix the value of d later) along 
A e and A' e , such that the two chains are separated by the supporting line of e. Note that the 
endpoints v and v' of e are the only ones that are not outside D e , and they are also in the 
kernel of D e (remember that A e and A' e can be chosen sufficiently flat). The edge centers are 
the only gadgets that have different edges in the source and in the target triangulation. Draw 
the edges that define the polygon Pu e in both T\ and T%. Let the first point of the upper 
chain be adjacent to all points of the lower chain in the source triangulation T\, and do the 
analogous for the last one in the target triangulation T2. We call the execution of flips that 
remove and add edges of an edge center transforming an edge center. 

3.1.2 Crossings 

If two straight-line edges e and / of G cross, also their corresponding circular arcs cross. Place 
one point at each of their four crossings (we will actually place them close to these crossings, 
see Appendix [AJ . In both source and target triangulation draw the edges connecting two 
points that are consecutive on any circular arc, which results in a crossing being represented 
by a convex quadrilateral, to which we add an arbitrary diagonal. 

3.1.3 Wirings 

Cf. Figure [9j Wirings are gadgets that represent the elements of V. Consider any vertex v 
of G and a small circle C with v in the embedding as its center. This point v is part of the 
triangulated point set. Place points on the crossings of C with the arcs of the edges incident 
to v in G. Since the graph is embedded on a convex n-gon and due to the small length of the 
vector t, these points occupy strictly less than half of C. This allows us to place two chains 
L and R, each of w points (the value of w to be defined later) on C in a way that any line 
between one point of R and one point of L separates v from the remaining construction. In 
both the source and target triangulation draw the edges between consecutive points on C. 
Draw a zig-zag path through the 2w points of L and R. Connect v to the two first points of 
the chains. The remaining part may be triangulated arbitrarily. 

The remaining faces in the two plane graphs we obtained so far are triangulated arbitrarily, 
however in a way that the resulting triangulations T\ and T2 have the same edges except at 
the edge centers. 

3.2 Analysis 

The basic idea of the construction is that an effective flipping algorithm has to choose which 
wirings to flip (requiring Aw — 2 flips each for flipping the edges of a wiring away and back 
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Figure 9: A wiring with its initial (and final) triangulation and a triangulation that allows to 
quickly transform the edge centers. 

again) in order that the triangulation of an edge center can be transformed using the point 
in its kernel at the chosen wiring. Also, the edges between and at the crossings need to be 
flipped away. We will fix the values of w and d to force this behavior of any flipping algorithm 
that requires less flips than a trivial upper bound. Every edge of G will be covered; using a 
vertex of G for covering corresponds to flipping all edges in the corresponding wiring. 

Let v and v' be any two adjacent vertices in G. The exact number of edges in T\ or T<i 
intersected by the segment vv' in the drawing might differ according to the choice of v and v' 
due to the triangulation of the wiring gadget at the region where the edge gadgets enter it. 
For every wiring there is at most a constant number c of such edges, giving in total at most 
t = cn edges. 

The following lemma shows how to deduce a flip sequence in our construction from a 
vertex cover of size k. Note that we do not claim that this is the optimum if k is optimal. 

Lemma 9. If there exists a vertex cover of size k in G, then there exists a flip sequence 
between T\ and T<i of length at most 

5 = 2(k(2w - 1) + m{Ax + 2d) + t) . 

Proof. Let C be a vertex cover of G with k = \C\. Let v be the vertex used to cover an edge. 
We use v to transform the edge centers of the adjacent edges in G (if they have not already 
been transformed). We need to flip all edges in the wiring to v, which takes 2w — 1 flips. Then 
we need at most c flips for the remaining wiring edges, as well as 2 further flips for the edges 
before the first crossing and 2 flips for the first crossing itself (all in all, with this method we 
need up to Ax + 2 flips for the crossing gadgets to "reach" the edge center). There, we need 
2d — 2 flips to make the edges incident to v (as in Figure [5]). Flipping in the desired way we 
would need at most 5 flips. □ 

On the other hand, a flip sequence should define a vertex cover. For the following lemma, 
we fix w > m(4x + 2d) + r + 1/2; further, we choose d such that (d — l) 2 > 2(n(2w — 1) + 
m(4x + 2d) + r) (note that since the term to the right is linear in d, such a value of d clearly 
exists and is polynomial in the problem size). 

Lemma 10. If there exists a flip sequence between T\ and T2 of length 5, then there exists a 
vertex cover of size k = , 5 2 . 
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Proof. We argue that the choice of d forces an effective algorithm to flip the wirings (which 
corresponds to covering vertices), and that the choice of w allows to transform the number of 
flips to the size of the corresponding vertex cover. 

If S > (d — l) 2 , then the choice of d implies a vertex cover of size re, which is trivially 
true. We therefore assume that 5 < (d — l) 2 . For every edge center, if we do not use the 
corresponding central points v or v' of a wiring, we would need at least (d — l) 2 flips due to 
Lemma [7j Suppose we want to transform an edge center D using a point v. As in the proof 
of Lemma [9j we need to flip all edges in the wiring to v, taking 2w — 1 flips. Note that this 
is optimal since only one of the edges can be removed with each flip. The values of d and w 
have been chosen in a way that flipping the edges of all wirings, crossings, and edge centers 
to the corresponding central point and back, as described, uses less flips than transforming 
one edge center, due to the bound of Lemma [9} For any algorithm, this means that nipping 
all edges at wirings and crossings twice and transforming the edge centers with a point at 
the wiring is cheaper than transforming one edge center without a point at a wiring. Due to 
Lemma [7] we know that we need a point at a wiring for each edge center to be transformed in 
fewer than (d — l) 2 flips. Therefore, we know that the (optimal) flip distance 5 op t is given by 

<>opt = &o P t(4w -2) + R for some R > . (1) 

Equation ([!]) shows how to deduce k opt from 5 opt ■ Note that if Aw — 2>R, the size of the 
minimum vertex cover can be calculated from the flip distance as stated. For the reasoning 
about approximation ratios later in this section, we actually want to have Aw — 2 > 2R. 
Lemma [9] gives us an upper bound on the flip distance, and hence R < 2(m(Ax + 2d) + r). 
Thus, we have chosen w such that nipping the edges of one wiring needs more flips than twice 
the number of all other flips in total in an optimal flip sequence. Therefore the flip distance 
5 opt gives 

^opt 

Aw - 2_ 

No matter how well an algorithm performs, it has to flip the wirings at least k op t times 
when using less than (d — l) 2 flips, and Lemma [9] tells us that the flips of the edges not in a 
wiring are in total fewer than those used for one wiring when the algorithm is optimal. □ 

The choice of w with respect to R and the proof of the previous lemma gives us the 
following fact. 

Corollary 11. For the flip distance 5 op t between T\ and T2 and the minimum vertex cover 
of size k opt we have k opt = \^_2 f or some positive R < 2w — 1. 

To show APX-hardness of the flip distance problem, we show that we have an AP- 
reduction [H pp. 256-261] from Cubic Minimum Vertex Cover using the previous lem- 
mata. Let k op t be the size of a minimum vertex cover for G and <5 op t be the flip distance 
between T\ and T2. The performance ratio of a minimization problem is the measure value of 
the approximation divided by the optimal value, e.g., k/k opt for an approximate vertex cover 
of size k. 

Definition 2 ([U pp. 257-258]). Let P\ and P2 be two NP optimization problems. P\ is 
AP-reducible to P2 if two functions f and g and a constant a > 1 exist such that: 



^opt 
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1. For any instance x of P\ and any rational r > 1, f(x,r) is an instance of P2. 



2. For any instance x of P\ and any rational r > 1, if there is a feasible solution of x, then 
there is a feasible solution of f(x,r). 

3. For any instance x of Pi and any rational r > 1, and for any y that is a feasible solution 
of f(x,r), g(x,y,r) is a feasible solution of x. 

4- f and g are computable by two algorithms whose running time is polynomial for any 
fixed rational r. 

5. For any instance x of P\ and any rational r > 1, and any feasible solution y for f(x, r), 
a performance ratio of at most r for y implies a performance ratio of at most 1 + a(r — 1) 
for g(x,y,r). 

In our case, / corresponds to the construction of the point set and the two triangulations. 
Requirements [T] and [2] follow from our construction. A vertex cover can be extracted from 
a flip sequence y from the edges flipped at the wirings; this analysis corresponds to g, and 
requirement [3] is therefore fulfilled. Both are polynomial-time algorithms (the parameter r is 
actually not used) , as demanded by requirement |4j 

Intuitively, Lemmata [9] and [10] give evidence that the reduction described so far fulfills 
also requirement [5] of Definition [2] However, there is a technical difficulty with the remainder 
term R that we have to take a closer look at. Let 5 be an approximate solution for the flip 
distance such that 5 < <5 pt r - Further, let R' be the remainder produced by the floor function 



in Lemma 10, i.e., 5 = i?'(mod Aw — 2). By Lemma 10, we get 



S-R' . S opt r - R' 



k < < 



4w -2 4w 



Let R be the remainder term for the optimal solution 5 op t as in Corollary 11 Then we have 

5 opt -R rR-R' , rR-R' , rR , r 

k < r — i = rk ont H < rfc ont H < rk (m t H — , (2) 

~ 4w - 2 Aw - 2 p iw - 2 - p Aw - 2 p 2 ' K ' 

where the equality and the latter inequality are due to Corollary |11[ Let a = 4 and suppose 
first that r — 1 = e > „, 1 , , . Then 

— 2fe pt + l 

r 1 e 1 / 1 \ 

rk opt + - = /c pt + aefeopt - 3efe op t + - + - = fc op t + aek opt + - - e I 3k opt - ~ 1 . (3) 

To get rid of the last part we use 

IV 3fc op t-l/2 1 



3kopt 2j- 2k opt + l > 2 
which, by ^ and ([3]), implies 

k < k opt + aek opt ■ 
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On the other hand, suppose that r — 1 = e < 2k 1 t+1 . Then from (2), we get 

r 1 e (I 

k < rk opt + - = fc op t + e£; opt + - + - = k opt + e I k opt + - 



r . , 1 e /, 1\ 1 

2 

< 



Ltll/2 1 
fe °pt + 2 fe opt + 1 + 2 " feopt + 1 ' 



Since the solutions to vertex cover are integers, this implies that k = k op t < k opt + aek op t. 
Hence, in both cases k < A; opt [1 + a(r — 1)] and our reduction fulfills all properties of an 
AP-reduction from Cubic Minimum Vertex Cover. 

Theorem 12. The problem of determining a shortest flip sequence between two triangulations 
of a point set is APX-hard. 



4 Conclusion 

In this paper, we showed that it is APX-hard to minimize the number of flips to transform 
two triangulations T% and T2 of a point set S into each other. As a by-product, Corollary [8] 
revealed an interesting aspect on distances in the flip graph. 

In presence of the recent NP-completeness result for simple polygons [2] , the main remain- 
ing open problem is the one for triangulations of convex point sets and its dual problem, the 
computation of the binary tree rotation distance [22j . 

Acknowledgements. The author wants to express his gratitude to Oswin Aichholzer, 
Thomas Hackl, and Pedro Ramos, as well as anonymous referees for valuable suggestions 
on improving the presentation of the result. 

A A Tedious Calculation of the Coordinates 

Section [3] already contained a description of the gadgets we used in our reduction. However, 
the validity of gadget-based reductions when proving NP- or APX-hardness for problems on 
point sets requires that the coordinates of the points used can be calculated in polynomial 
time. While this is often, implicitly or explicitly, left to the reader, we perform this step in 
high detail in this appendixj^] 

The reader might have noticed that our high-level construction involves points placed at 
the crossing of circular arcs, which, in general, leads to irrational coordinates, even if the 
circular arcs are defined by rational points. We will give a construction that slightly varies 
from the one described that uses only rational coordinates, with both the numerator and 
denominator bounded by a polynomial in the input size. 

The gadgets in our reduction do not make use of collinear points. However, we did not 
explicitly mention how to avoid three points on a line when the construction is assembled. In 

2 This was explicitly demanded in a review of a conference submission of this paper. 
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Figure 10: Five points on the unit circle; a point at the red cross would introduce a supporting 
line through a crossing of two other supporting lines and is therefore forbidden. 



this appendix we give an explicit construction of the point set in a way that it is in general 
position, i.e., that no three points are collinear. 

Note that the construction might not be "economic" in the sense that the construction 
might be possible with smaller values. We will always prefer constructions that are easy to 
prove. We will place the points on and close to the unit disc (meaning that a coordinate will 
never exceed 2), and therefore can give the size of a coordinate in terms of the size of its 
denominator. 



A.l Placing the Points of the Convex Polygon 

As a first step, we give a simple construction of a convex n-gon for placing the central points 
of the wiring gadgets with all vertex coordinates being rational and the denominators being 
in 0(n 10 ). Further, we want to assure that no three diagonals cross in the same point. For 
doing so, we will first choose n 5 candidate points on the unit circle and then select n points 
out of them. 

Rational points on the unit circle are known to be given by ^fzq^i; with t E Q, 

see, e.g., [TO]. We define a sequence K of candidate points with t = i/n 5 for the integers 
1 < i < n 5 . (For consistency with later parts and ease of presentation therein we choose the 
candidate points from the upper-right quadrant in counterclockwise ascending order; hence, 
the value of t is between and 1.) Now we select n points out of K such that there are no 
three diagonals that cross at a single point. We choose the first five points of our final set 
from the candidate points. Suppose we have chosen j > 5 points such that no three diagonals 
cross at a single point. We have n 5 — j points in K to choose the next point from. Look at all 
possible combinations of five points among the already chosen ones. Each combination gives 
exactly five points on the unit circle that can not be chosen, and this point can be none of the 
j already chosen candidate points. Hence, we have 5(g) + j "forbidden" points (which might 



not be among the candidate points). See Figure 10 However, we have n 5 > > 5(g) + j 
candidate points to choose from, and therefore we for sure can choose point number (j + 1). 
We denote this set of points by Py; the elements of Py are the points representing the vertices 
of the input graph. 
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IV 



Figure 11: Construction to obtain bounds for 5. 



Proposition 13. A point set of n points in convex position with all coordinates rational 
having a denominator in 0(n 10 ) and no three diagonals crossing in the same point can be 
found in polynomial time. 

Note that the facts that no three diagonals of the resulting n-gon intersect and that the 
coordinates are bounded also give us a lower bound on the distance between intersection 
points and other diagonals, which we will use in the next part. 



A. 2 A Sufficiently Small Value 

Find the minimum squared distance from every of the Q) crossings of the diagonals of the n- 
gon to the diagonals not involved in the corresponding crossing. Since this squared distance 5 2 
is given by (x a — Xb) 2 + (y a — Ub) 2 between two points a and 6, we can set 5' e = \x a — x&| to 
obtain a "small", rational and positive distance 5' e < 5 e (at least one of the horizontal or 
vertical distance is non-zero, in particular, up to here no two points can have the same x- or 
y-coordinate) . When we construct the tunnels that are formed around an edge of the drawing 
of the input graph, we can choose, say, S' e /3 as an upper bound for the distance between the 
edge and the edges defining the tunnel. Then no three tunnels have a non-empty intersection. 
(Our actual tunnels will be even narrower.) 

The vertex gadgets used "small" circles around each point in Py. Let u, v, w be a consecu- 
tive triple of vertices on the n-gon defined by Py . Let 5 2 denote the smallest squared distance 
between v and the line through u and w for every choice of the triple. As with the tunnels, we 
can choose a rational 6' v < 5 V by choosing only the horizontal or vertical distance between v 
and the closest point on the supporting line of u and w. Let n ov be the line through v that is 
perpendicular to the supporting line of the origin o and v. Consider the distances from u and 
w to n ov . Let 5 2 be the smallest squared distance for all choices of v (and corresponding u 
and w), and choose a rational 8' < 5 n as before. Further, let 5 r be the smallest horizontal or 



vertical distance between two points in Py (which is non-zero by construction). See Figure 11 
We define 5 = min{5' e /3, 6' v , S' n , 5 r }. If we now choose the radius of the cycle centered at each 
vertex by ry = 5/6, then no two circles intersect (there is actually a distance of at least 4ry 
between two circles), and each circle only intersects the edges of the input graph that are 
incident to the vertex it is centered at. Further, no circle intersects the convex hull of two 
other circles. 
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A. 3 Tunnel Construction 



For each edge e of the input graph connecting two vertices v and w, we now give the construc- 
tion of the tunnels. Let C v and C w be the circles around v and w, respectively. The tunnel 
for the edge between v and w is given by two segments, each having one endpoint on C v and 
one endpoint on C w . We want to get rational points on C v and C w . Since these cycles are 
not only defined by a rational center point, but also have a rational radius, the problem boils 
down to find a rational point on the unit circle, or, equivalently, a (possibly irrational) angle a 
such that sin(a) and cos(a) are rational, within some interval given by quadratic irrationals. 
Sines with this property are called rational sines, and correspond with the parametrization 
of the unit circle that we already used before. Canny, Donald and Ressler [10] give an algo- 
rithm for finding a rational sine for a parameter t = p/q, such that \p/q — x\ < e, for given x 
and e (we will use an extended method for non-rational radii later). Their algorithm gives a 
denominator q in 0(l/e), and the running time is polynomial in q. However, the input x is 
an approximation as well, and their goal is to get rational sines with small binary represen- 
tation. Our angle intervals, however, are given by rational points and their relative position 
to the circle center. For finding a point within this interval, the Farey approximation as used 
in [TO] for t = p/q is sufficient and easy to apply for our setting, as we do not need an explicit 
approximation of the angle and the interval as input (this algorithm searches a point inside 
the interval in the fashion of binary search, calculating the mediant of two rational values in 
each step). We, however, need an upper bound on the denominator q derived from the points 
defining the angle. 

Now we show how to use the results of [10J for our needs. Consider the unit circle and 
two points a and b. Let Za and Zb be the polar angles of these points, and w.l.o.g. let 
Za < Zb. We describe only the case where both angles are within [0, . . . , tt/2], the other cases 
are similar (and can easily be distinguished). To approximate an angle between Za and Zb 
using a rational number t, we reason about the (possibly irrational) values t a and t&. For t a 
and Za we define 

It 

• i / \ a 
sm(Za) = ^ , 

which, when choosing the appropriate root, gives 



sin(Za) y sin 2 (Za) 



The sine of Za is given by a y j \Ja% + a*. We define the analogous for tj. The values of Zb and 
tb are defined analogously. We therefore need to find a rational number t with t a < t < £&. The 
Cauchy bound (see [24J ) for an algebraic number g being the root of a polynomial YliLo c ^ x% 
with rational coefficients Cj is given by 

i i ^ l c °l 
\9\ > 



c | + max{|ci|, . . . , \c m \} 



The difference \t a — tb\ is therefore bounded from below by a rational that has a denominator 
polynomial in the problem size. Using Farey approximation, we can find a rational t whose 
denominator exceeds the denominator of the bound only by a polynomial factor. 
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Figure 12: Construction for the tunnel endpoints: The two extremal candidate points for p v 
are chosen inside the two gray wedges. Note that v and w in the drawing do not fulfill the 
required vertical distance since the drawing would get too small. 



Since we can choose rational points on the unit circle inside an interval (and therefore on 
instances of the unit circle that are translated and scaled by rational values), we now have the 
tools to choose the endpoints of the tunnels. For two vertices v and w, let these be called p v 
and q v (placed on C v ), as well as p w and q w (placed on C w ). Hence, a tunnel between v and w 
consists of the quadrilateral p v q v q w p w . In order to prevent collinear triples of points, we again 
select a set of candidate points on C v and C w and choose the four points among them. Note 
that this results in tunnels that might not be exactly rectangular, but this is irrelevant for 
our final construction. See Figure 12 for an accompanying illustration. W.l.o.g. suppose that 
x v < x w . To obtain the set of candidate points for p v , consider the segment between w and 
the point (x w ,y w + ry), where ry is the radius of the circles, which we call the upper spoke 
of w. Let the lower spoke of w be defined analogously between v and the point (x w ,y w — ry). 
Find the two parameters t\ and t<i for rational points pt ± and pt 2 on C v such that the line 
through v and pt x intersects the upper spoke of w at a point in the topmost eighth of the 
upper spoke and the line through v and pt 2 intersects the upper spoke in its fifth eighth from 
below. We can now select our set K v of candidate points from the interval [£i, ^2] • The same 
can be done for two parameters t% and £4, with the roles of v and w changed. We select a point 
p v £ K v and a point p w 6 K w as the endpoints of one side of the tunnel gadget between v and 
w. Note that we do not need to require the sides to be parallel to the supporting line of v and 
w (we could do so by increasing the number of candidate points). The crucial property of the 
points we need is that p v vwp w forms a convex quadrilateral and we therefore have to prove 
that p v is always left of the directed line through v and p w (and, analogously, that p w is right 
of the directed line through w and p v ). Let l w = (x w , y w + ry/2) and u v = (x v ,y v + ry). The 
diagonals vl w and u v w of the trapezoid u v vwl w intersect each other at a ratio of (ry/2)/ry, 
i.e., at two thirds of the interval [a^,,a; w ]. Recall that the radius ry was chosen in a way that 
the disc centers have a horizontal distance of at least 6ry. Hence, the segments intersect 
outside C w ; the topmost candidate point on C w is below the line through v and the lowest 
candidate point on C v , and vice versa. 
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Note that since the candidate points on C v are chosen inside the convex hull of C w and 
v, no two tunnels from v can intersect. 

It remains to find the correct number of candidate points. Suppose we already constructed 
all but one tunnel point. Since at every circle there are at most six tunnel points there are less 
than ( lines on which we are not allowed to place a point. Every line intersects the circle 
we place the last point on at most twice. Hence, if we choose more than twice the number of 
points as we have lines, we can always choose a point to remain in general position. 



A. 4 Points in the Tunnels 

For each tunnel, we construct two circular arcs, one for each segment defining the tunnel, 
on which we place the points of the edge center. The crucial property of such an arc is that 
for two wire centers v and w, these points are the only ones in the kernel. Let q v and q w be 
the two endpoints of a tunnel edge, s.t. q w is to the right of q v and the interior of the tunnel 



is above the line q v q w . See Figure [13) The constructed arc will start at q w and end at q v . 
W.l.o.g. let the angle between q v q w and q v w be smaller than the one between q v q w and q w v. 
Then we consider four rays, namely the ones that leave q v to the right in an angle of 0, (— 7r/4, 
and 7r/4 with the x-axis and the upward vertical ray at q v . Let r be the one that opens the 
smallest positive angle a with q v q w . If the angle between q v and q w and q v w is smaller than 
a, then let s be the ray through w starting at q v ; otherwise, let s = r. Construct the circle A 
that has its center on the line through the midpoint of q v q w perpendicular to q v q w and that 
has s as a tangent. The coordinates of the center of A are still rational. It is well-known that, 
when given any rational point p on A and a line t with rational slope that intersects A at p 
and a second point p', the point p' is rational as well, see, e.g., [Ml p. 5]. Hence, we need to 
appropriately choose lines through a point p. The crossings of the segments that define all 
the tunnels identify the region where to place the edge center. Let R be the region we have 



to place the points in (marked gray in Figure 13). By the choice of r, we constructed A in 



a way that we can mirror and rotate the plane orthogonally such that the intersection of A 
and R is within an angle of and 7r/4 from q v . This means that any line £ through q v and 
this intersection will have a slope t between and 1. This reasoning is similar to the one of 
Burnikel [9] to adapt the techniques of [10] for such rational circles (i.e., circles given by three 
rational points). As before, we can use, e.g., Farey approximation for the slope t of £. At 
each iteration, we check whether the second intersection of I with A is inside the quadrilateral 
R, and, if not, on which side it is. Since the denominators of the coordinates of the points 
defining A and R are polynomial, there is a polynomial lower bound on the difference between 
the (possibly non-rational) parameters for the two points where A enters and leaves R (as for 
the construction of the tunnel endpoints). Hence, after a polynomial number of steps, we have 
a rational slope for I such that I passes through A inside R; therefore, also this intersection 
point has rational coordinates and its denominator is polynomial in the problem size. To 
obtain a second such point, the process can be continued. Now we have two points in the 
intersection of A and R which define two slopes of lines through q v . Any line through q v with 
a slope in the interval between these two slopes gives a rational point on A n R. Hence, we 
can choose our candidate points by dividing that interval. 

The remaining problem is the one of choosing the points for the crossing gadgets. Two arcs 
in crossing tunnels will, in general, cross at a point that does not have rational coordinates. 
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Figure 13: We want to choose rational points on the (blue) arc inside the gray region by 
Farey approximation on the slope t. Note that the gray region actually is, for presentational 
reasons, drawn too close to w. 




w 



Figure 14: Construction of tunnel crossings. The drawing shows the lower part of a tunnel 
between v and w and a part of another tunnel (indicated by the near- horizontal strokes). 



The crucial property of the points of the crossing gadgets, however, is that they are outside the 
edge centers (recall Definition [l]) and that the edges between them can "quickly" be flipped 
to the center of the corresponding wiring gadget. Placing the points for the crossing gadgets 
at the crossings of the segments that define the tunnels would satisfy these constraints, but 
would lead to collinear triples. So we have to slightly perturb each point p to obtain a point p' 
without loosing these properties. See Figure 14 Between every consecutive pair of crossing 
points on a tunnel segment QvQw we can choose the rational midpoint. If the perturbed point p 
remains on the same side of the line through the wire center and the midpoint as p, the order 
around the wire center is maintained. Further, the perturbed points have to remain on the 
same sides of the lines that define the wedges of the edge centers involved. Together with the 
tunnel edges, these constraints give a convex region from which we can choose our perturbed 
point. We may again place a circular arc inside this region (marked gray in Figure 14) on 
which we select a sufficiently large number of candidate points, analogously to the construction 
of the other gadgets. 



A. 5 Points for the Wiring 

Finally, we place the points at the wiring gadgets that allows to draw the wiring edges, see 
Figure [15] The circles for the wiring gadgets are scaled versions of the unit circle. For a 
vertex v, let n ov be the line through v that is perpendicular to the supporting line of the 
origin o and v. Since the coordinates of v are rational sines, the intersection points of n ov 
with the circle C v are rational as well. Due to the choice of S' n , all points on C v that define 
tunnels are on the same side of We are given two intervals, each between two rational 

sines, i.e., between the "extremal" tunnel endpoints on C v and the intersection points of n oe 
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Figure 15: Construction of the wiring points. The small gap (indicated by the arrow) on the 
circle C v of v between the intersection point with n ov and the neighboring tunnel endpoint 
can be used for the candidate points. 

with C v . Therefore, we can choose a sufficient number of rational candidate points on C v to 
choose the points for the wiring from. 

A. 6 Concluding Remarks on the Embedding 

The crucial part throughout the whole embedding procedure is that each (intermediate) point 
that is not a candidate point is constructed using only a constant number of other points. 
The candidate points were constructed with polynomial parameters. Hence, all denominators 
are polynomial in the input. In particular, note that even though some intervals were defined 
by points with algebraic coordinates, a lower bound on the interval can be given in terms of 
the other, rational coordinates that were used in the construction. This allowed us to find 
rational points with polynomial denominators within these intervals. 
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